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TRANSONIC  WING  REDESIGN  USING 
A GENERALIZED  FICTITIOUS  GAS  METHOD 


I . SUMMARY 


A numerical  method  for  transonic  shock-free  or  nearly  shock-free  airfoil 
and  wing  redesign  based  on  the  full  potential  equation  is  presented.  The 
method  utilizes  a generalized  fictitious  gas  approach  wherein  a variety  of 
parameters  controlling  the  character  of  the  fictitious  gas  laws  are  introduced 
to  provide  a degree  of  control  over  the  redesigned  upper  surface  geometry  and 
the  pressure  distribution  of  the  redesigned  shape.  Results  for  a redesigned 
advanced  airfoil  as  well  as  a three-dimensional  wing  are  illustrated. 
Significantly  improved  aerodynamic  characteristics  are  achieved  through  the 
present  redesign  technique. 


II.  INTRODUCTION 


It  is  well  known  that  the  subsonic  cruise  speeds  of  high  performance 
aircraft  are  limited  by  the  onset  of  the  transonic  drag  rise.  Such  drag  rise 
is  primarily  caused  by  the  formation  of  shock  waves  which  recompresses  the 
local  supersonic  flow  to  subsonic  flow.  One  means  for  delaying  the  drag  rise 
and  its  associated  unfavorable  flow  phenomena  is  to  improve  (redesign)  the 
wing  contour  shape  so  as  to  delay  the  formation  of  strong  shock  waves. 

Several  transonic  wing  design  concepts  have  been  explored  over  the  years. 
The  peaky  airfoil  concept  first  introduced  by  Pearceyl  leads  to  the  design 


of  stable  shock-free  transonic  airfoils.  However  this  peaky  type  of  airfoil 
usually  does  not  display  acceptable  off-design  characteristics.  Experimental 
results  indicate  that  at  lower  Mach  number  the  high  pressure  peak  near  the 
leading  edge  can  trigger  early  boundary  layer  separation  and  cause  high  drag 
levels.  On  the  other  hand,  a design  with  a roof-top  type  pressure 
distribution  may  be  favorable  to  the  boundary  layer  flow,  but  it  may  not  be  a 
very  stable  shock-free  design. 

Current  design  efforts  include  the  hodograph-based  works  of  Nieuwland?, 
Bauer,  Garabedian,  and  Korn^,  Boerstoel4  and  Sobieczky^,  and  the  use  of 
optimization  techniques  (Hicks,  Murman,  and  Vanderplaats6)  for  the  design  of 
supercritical  airfoils  without  drag  creep.  Iterative  finite  difference 
methods  have  been  explored  by  Steger  and  Klineberg^,  Tranen®,  Carlson^, 
Shankar,  Malmuth,  and  Cole^H^  and  Volpe^,  wherein  the  design  problem 
is  treated  as  an  inverse  problem,  i.e.,  the  airfoil  shape  that  supports  a 
given  pressure  distribution  is  found  by  iteration.  More  recently,  Sobieczky 
et  al^#>  an(j  Yu*4  have  developed  a unique  fictitious  gas  method  for  the 
design  of  shock-free  airfoils  and  wings. 

The  present  research  extends  the  previous  works^.14  to  a general 
airfoil  and  wing  redesign  method.  The  method  utilizes  a generalized 
fictitious  gas  approach  in  the  redesign  process  wherein  a variety  of 
parameters  controlling  the  character  of  the  fictitious  gas  laws  are  introduced 
to  provide  a degree  of  control  over  the  redesigned  upper  surface  geometry  and 
the  pressure  distribution  of  the  redesigned  shape.  The  gas  dependent 
redesigned  configurations  can  be  systematically  varied  to  explore  the  trades 
between  aerodynamic  performance  and  geometry  (thickness)  changes  in  the 
transonic  regime  of  flight. 
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III.  FORMULATION  OF  GOVERNING  EQUATIONS 


The  basic  equations  to  be  used  in  the  present  redesign  procedure  are  the 
conservation  of  mass,  energy  and  entropy  for  a potential  flow,  i.e., 


(P<f>x)x  + + C/V* )j,  =o. 


j _ _L  4 J'-L. 

a t i 2 , 


f - ( C a) 


2 \ ~ 1 


respectively,  where  <$>  is  the  complete  velocity  potential, P,  a,  and  q are 
the  conventional  definitions  of  the  local  density,  speed  of  sound  and  flow 
speed  respectively.  All  flow  variables  are  normalized  with  respect  to 
freestream  conditions. 


The  present  redesign  technique  requires  modification  of  a ''“liable,  fully 
conservative  analysis  code  for  the  calculation  of  transonic  flow  over  an 
initial  input  configuration.  In  two-dimensional  redesign,  one  can  either  use 
Jameson's  full  potential  code^  or  Holst's  approximate  factorization  (AF) 
code^  as  a basis;  where  in  three-dimensional  redesign,  we  use  a modified 
version  of  Jameson-Caughey's  finite  volume  code^  for  our  redesign  method. 


L 
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IV.  REOESIGN  PROCEDURE 


The  present  redesign  method  begins  with  the  geometry  of  an  initial  input 
configuration  that  typically  produces  a shock  wave  terminating  the  upper 
surface  supersonic  bubble.  The  method  then  provides  a local  modification  of 
the  upper  surface  geometry  underneath  the  supersonic  bubble  so  as  to  produce  a 
shock-free  or  nearly  shock-free  flow. 


The  principal  redesign  concept  embodies  a temporary  modification  of  the 
governing  equations  locally  in  the  supersonic  region  so  as  to  retain  an 
elliptic  behavior  within  the  supersonic  bubble-hence,  the  terminology 
"fictitious  gas".  In  two  dimensional  flow,  it  can  easily  be  proved  that  the 
following  density  laws  will  ensure  elliptic  behavior  of  the  governing 
equations  in  the  supersonic  region.  One,  termed  a continuous  density  law,  is 


with  both  and  Cj  greater  than  zero.  Another  is  termed  a discontinuous 
density  law,  meaning  that  the  first  derivative  dP/dq  is  discontinuous  at  the 
sonic  surface,  and  is  given  by 


Pz 


(5) 

U 


with  > -1,  C.?  ^ 0.  In  both  equations  the  superscript  refers  to 
the  sonic  condition.  Figure  1 illustrates  the  density  variations  versus  q/a* 
for  several  Pj,  Cj  and  P?,  C?. 
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In  three-dimensional  flow,  only  the  discontinuous  density  law,  equation 
(5),  with  proper  combination  of  Pp  and  C£  has  been  used.  More  general 
density  laws  remain  to  be  explored. 

The  governing  equations  and  the  boundary  conditions  in  the  subsonic 
regions  remain  unchanged.  Numerical  solutions  to  these  modified,  elliptic 
type  governing  equations  contain  no  shock-wave-like  discontinuities,  and  are 
correct  in  the  subsonic  regions  up  to  the  sonic  boundaries  that  separate  the 
subsonic  and  supersonic  regions.  The  results  within  the  supersonic  bubble 
where  the  gas  is  "fictitious"  are  discarded. 

The  results  along  the  sonic  surface  produced  by  the  fictitious  gas 
calculations  are  used  as  initial  data  for  the  recalculation  of  the  flow  field 
in  the  supersonic  region  using  the  real  density  law.  In  the  two-dimensional 
airfoil  redesign  procedure,  both  a characteristics  method  and  a marching 
method  have  been  used,  while  in  the  three-dimensional  study,  only  the  marching 
method  is  used  to  recompute  the  continuous  supersonic  flow. 

The  characteristics  method  used  in  the  two-dimensional  supersonic  flow 
calculation  has  been  detailed  in  reference  13  and  will  only  be  briefly 
discussed  here.  Basically,  the  calculation  is  done  in  a hodograph-1 ike 
working  plane,  in  which  the  characteristics  are  orthogonal  straight  lines. 

The  velocity  potential  and  the  stream  function  Ip  are  solved  along  the 
characteristics  net,  while  the  physical  coordinates  of  each  characteristic 
line  are  obtained  from  a simple  integral  relationship.  Finally,  the  new 
redesigned  airfoil  geometry  underneath  the  supersonic  bubble  that  provides  the 
continuous  supersonic  flow  is  obtained  through  interpolation  of  zero  stream 
function  in  the  flow  field.  Because  the  mass  flux  is  exactly  conserved  in  the 
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difference  aoprox imation,  the  redesigned  geometry  that  produces  a shock-free 
solution  is  smooth  and  continuous  to  the  accuracy  of  the  finite  difference 
approx imations. 

The  marching  procedure  which  is  used  both  in  the  two-dimensional  airfoil 
and  in  the  three-dimensional  wing  redesign  is,  in  general,  less  accurate  in 
numerical  resolution,  but  is  more  flexible  and  easier  to  implement.  Because 
the  sonic  surface  data  generated  from  the  fictitious  gas  analysis  are  smooth 
and  continuous,  one  can  obtain  the  first  derivatives  of  the  flow  variables  on 
the  surfaces,  i.e.,  Ux,  Vx,  Wx,  Uy,  V , W and  Uz,  Vz,  Wz  through  the 
use  of  Taylor  series  expansion  of  U,  V,  W along  two  arbitrary  directions  on  the  sonic 
surface,  plus  the  continuity  equation  in  quasilinear  form 

(a2-u2)^  -t(cL*  + (.&*-&*)  <ujz 

-2  UTS  UH  ~ 2 VOf  tfz  -2UiOU£  = 0, 

(6) 

and  the  three  irrotationalit.y  conditions 


ut 

= Vy9 

(7a) 

^ l 

(7b) 

■=  rfy. 

(7c) 

Once  the  derivative  terms  are  calculated,  one  can  then  compute  U,  V and  W on  the  next 
layer  of  computational  surface  through  a direct  numerical  integration  procedure. 

Here,  as  an  illustrative  example,  we  choose  the  coordinate  Z - the  direction  normal  to 


- 


. . . .. 
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the  wing  surface  - as  our  marching  direction.  The  flow  variable  on  the  next  layer 
of  computational  surface  is  obtained  by  the  trapezoidal  integration  formula 

U2  ^ U,  + 0.5  ( -+  Hg2)i  (8) 

here,  subscript  "1"  and  "2"  refer  to  the  old  and  new  computational  surfaces 
respectively.  The  marching  procedure  continues  until  the  whole  supersonic  region  is 
recomputed.  The  new  wing  shape  and  the  pressure  distribution  associated  with  it  are 
then  extracted  from  the  flow  field  solution.  Here,  the  stream  surface  condition 

U,  -w  -f  i rZy  ^ o,  (9) 

is  used  to  iterate  the  new  wing  surface  coordinates.  A more  detailed  formulation  for 
this  marching  technique  can  be  found  in  reference  14. 

The  final  step  of  our  redesign  procedure  consists  in  verifying  the  redesigned 
configuration  to  ensure  that  it  does  indeed  provide  a shock-free  or  a nearly 
shock-free  pres:ure  distribution.  In  two-dimensional  airfoil  redesign,  the  off-design 
charcteristics  of  the  resultant  airfoils  have  also  been  studied. 


V.  NUMERICAL  RESULTS 


Two-Dimensional  Airfoil  Redesign 

An  advanced  wing  section  at  freestream  conditions  * .75,  ci  - 1° 
is  selected  as  an  illustrative  example  for  the  present  redesign  technique. 
Figures  2a-c  show  the  sonic  line  shapes,  the  surface  geometries  and  the 
pressure  distributions  of  the  original  (curve  0)  and  the  redesigned  (curves  1 
to  5)  airfoils.  Here  curves  1 to  5 correspond  to  the  redesigned  results  of 
the  following  density  formulations  in  the  supersonic  region: 

Table  1.  Density  laws  used  in  the  redesign  procedure 


Curve 

density  equation 

P1  or  P2 

Cj  or  C2 

1 

4 

1 

0.5 

2 

4 

8 

0.5 

3 

5 

-0.5 

0 

4 

5 

0.0 

0 

5 

5 

1.0 

0 

The  density  laws  employed  in  the  fictitious  gas  analysis  exlblt  great 
control  over  the  shape  and  the  size  of  the  sonic  surface,  which  In  turn, 
provide  a variety  of  choices  for  the  redesigned  airfoil  shape  as  well  as  Its 
associated  pressure  distribution. 
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The  off-design  characteristics  for  the  redesigned  airfoils  are  plotted  in 
Figure  2d,  where  the  lift  and  drag  coefficients  are  calculated  at  a fixed  Mach 
number.  Notice  that  there  is  a systematic  correlation  between  the  decreasing 
thickness  of  the  airfoils  and  the  inviscid  drag  improvements.  Figures  2e  and 
2f  illustrate  the  off-design  pressure  distributions  for  the  original  airfoil 
and  the  redesigned  airfoil  case  3.  The  redesigned  shape  produces  a pressure 
distribution  with  weaker  shock,  and  consequently,  lower  wave  drag  value  for  a 
range  of  lift  coefficients.  The  effect  of  the  redesigned  pressure  field  on 
the  boundary  layer  flow  has  not  been  investigated  although  this  can  easily  be 
done  by  a simple  interaction  with  a boundary  layer  code.  For  practical 
airfoil  redesign,  such  interaction  calculations  would  be  important. 

The  present  redesign  procedure  can  also  be  readily  implemented  in  other 
fully  conservative  analysis  codes.  To  demonstrate  that,  we  modified  Holst's 
approximate  factorization  (AF)  code  for  the  fictitious  gas  analysis  and 
repeated  the  redesign  of  the  same  airfoil  using  the  constant  density  law 
in  the  supersonic  region.  The  sonic  surface  shape  and  the  redesigned  airfoil 
geometry  obtained  from  this  code  are  compared  with  that  of  Jameson's  full 
potential  code  in  Figures  3a  and  3b.  The  two  results  are  practically 
identical . 

The  present  two-dimensional  redesign  procedure  is  extremely  efficient. 

The  complete  redesign  process  requires  only  one  third  of  the  CPU  time  of  a 
corresponding  transonic  flow  analysis  calculation.  This  is  not  suprising, 
since  the  elliptic  type  of  governing  equation  can  be  effectively  solved  by  a 
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Three-Dimensional  Wing  Redesign 

The  three  dimensional  redesign  results  illustrated  here  provide  a 
demonstration  of  feasibility  of  the  present  redesign  technique.  Figures  4a-c 
illustrate  the  original  and  redesigned  results  for  a nonlifting  aspect  ratio 
1,  rectangular  wing  with  an  NACA  64A  series  airfoil  at  = .925.  The 

original  thickness  ratio  tapered  linearly  from  GX  chord  at  the  root  to  .6% 
chord  at  the  tip  so  as  to  generate  a closed  sonic  surface  inboard  of  the  tip. 
The  original  configuration  showed  a rather  strong  shock  at  the  inboard 
station,  an  expected  phenomenon  for  a conventional  NACA  6-digit  wing  section 
at  a high  subsonic  speed.  The  redesigned  results  using  the  density  laws  3,  4 
and  5 in  Table  1 are  plotted  on  the  same  figures.  Again,  the  use  of  different 
density  laws  in  the  redesign  procedure  effectively  controls  the  redesigned 
shapes  and  their  associated  pressure  distributions,  although  their  effects  are 
slightly  less  pronounced  than  in  the  two-dimensional  studies.  Both  the  shock 
wave  and  the  wave  drag  have  been  effectively  eliminated  by  the  present 
redesign  procedure. 

The  computational  efficiency  for  the  three-dimensional  redesign  procedure 
is  less  effective  than  that  of  the  two-dimensional  redesign  procedure.  This 
is  simply  because  we  have  not  utilized  the  elliptic  property  of  the  nodified 
governing  equation.  Currently,  the  three-dimensional  redesign  calculation 
requires  about  the  same  amount  of  CPU  time  as  that  of  a typical  transonic 
analysis  calculation.  Improvement  on  the  computational  efficiency  utilizing 

1 O 

accelerating  schemes,  such  as  Holst  and  Ballhaus'  AF  method  , would  be 
desirable  for  the  further  exploration  of  the  present  redesign  method. 


The  present  redesign  method  does  not  guarantee  the  existence  of  a shock- 
free  solution  for  an  arbitrary  initial  input  configuration.  In  fact,  both  in 
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two-dimensional  and  in  three-dimensional  redesign  calculations,  limit  lines 
may  occur  for  cases  with  a large  and  steep  supersonic  region.  In  such  a case, 
the  initial  data  along  the  sonic  surface  may  be  overspecified,  and  a 
continuous  shock-free  flow  could  not  be  obtained.  Practically,  this  does  not 
cause  any  serious  problem,  since  an  approximate  solution  which  produces  weak 
shock  can  usually  be  constructed  from  the  present  redesign  method. 


VI.  CONCLUSION 

A unique  and  practical  transonic  shock-free  or  nearly  shock-free 
configuration  redesign  method  has  been  explored.  The  method  utilises 
generalized  fictitious  gas  laws  to  control  the  shape  and  the  size  of  the 
supersonic  region  of  a given  initial  configuration.  The  gas-dependent  sonic 
surface  determines  the  character  of  the  redesigned  configuration.  The  present 
method  is  especially  effective  for  airfoil  or  wing  redesign  which  is 
restricted  to  a minor  modification  of  the  upper  surface  geometry  to  achieve  a 
favorable  transonic  pressure  field.  The  method  might  also  be  particularly 
advantageous  for  variable  geometry  (adaptive)  wing  design  work  where  it  is 
desirable  to  obtain  shock-free  or  nearly  shock-free  flow  over  a range  of 
flight  conditions  by  variation  of  the  wing  surface  geometry.  With  some 
additional  programming  effort,  the  method  can  be  extended  to  the  redesign  of 
wing-body  combinations. 
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Real  density  (curve  0)  versus  fictitious  density 
(curves  1 to  5) 


Figure  2e.  Off-design  pressure  distributions  for  the  original 
§ airfoil  at  M*  = .75 


Figure  2f. 


Off-design  pressure  distributions  for  the  redesigned 
airfoil-case  3 it  M.  * .75 
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Figure  3a.  Sonic  surfaces  of  an  advanced  wing  section  at  = .75 
and  oi-  1°  generated  by  the  modified  Jameson's  and 
Holst's  codes  using  constant  density  law  - curve  4. 
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Figure  4c.  Original  (curve  0) 
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